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INIROBUCIION 


The alternation of cold and warm stages of glacial and inter- 
glacial periods during the Pleistocene had a strong impact on the 
distribution and evolution of vegetation and associated faunal spe- 
cies (COMES & KADEREIT 1998; Hewitt 2000). One of the animal 
genera strongly influenced by the succession of glacial and intergla- 
cial stages has been the plethodontid salamander Speleomantes, due 
to its particular biological features. Indeed these salamanders usually 
hatch completely metamorphosed from eggs, lack an aquatic larval 
stage, are lungless with cutaneous and bucco-pharyngeal respiration 
and are strongly associated with underground retreats, offering a 
steady cold and moist environment (SALVIDIO et alzi 1994; CIMMA- 
RUTA et alu 1999; Lanza 1999). These features, coupled with low 
vagility and restricted home ranges, make Speleomantes’ geographic 
distribution strongly affected by climatic variations. 


In particular, the species Speleomantes strinatit has a relatively 
wide range, compared with those of the congeneric species, com- 
prising north-western Italy and south-eastern France. In these zones 
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the climatic variations of the Quaternary had great influence on the 
evolutionary history of local populations because of the uneven 
local orography. The genetic diversity of this species, investigated 
by means of allozyme electrophoresis, has been correlated with the 
paleoclimatic events of the species range, to understand the mi- 
croevolutionary mechanisms leading to its extant genetic structure. 
A number of subspecies has been described on morphological basis 
within this species (see LANZA et alu 1995: p.53) but their validity 
has been recently challenged based on genetic ground (NASCETTI et 
alti 1996). The data obtained here on the genetic—geographic struc- 
ture of S. strinati have been used to reanalyse the systematics of 
these questioned taxa. 7 


MATERIALS AND METHODS 


About 400 specimens of S. strinati, collected from 29 popula- 
tions scattered through the whole range of the species have been 
studied (Tab. 1). 


The techniques used for multilocus electrophoresis are those 
reported in LANZA et ala (2005) and are described in detail in 
LANZA et ala (1995). Following these procedures with minor modi- 
fications, the following 33 putative loci have been scored: a—-Gpdh, 
Ldh-1, Ldh-2, Hbdh, Mdh-1, Mdh-2, Mdhp-1, Mdhp-2, Idh-1, 
Idh-2, 6-Pgdh, Gapdh, NADH-dh, Sod-1, Sod-2, Np, Aat-1, 
Aiat=2 ClgeAclie,Bigt=4,vPep=C, Ped-D, Lap, Ada-1,.Ada-2, Ca-2, 
Ca-3, Fum, Mpi, Gpi, Pgm-1, Pgm-2. They included about one 
third of fast and two thirds of slow evolving loci, sensu SARICH 


(1977) and KoEHN & Eanes (1978). 


The genetic variability of populations was evaluated using the 
following parameters: overall number of alleles per population (4); 
percentage of polymorphic loci (P) following 99% (P99) and 95% 
(P95) criteria and mean expected (He) and observed (Ho) heterozy- 
gosity per population. The statistical significance of departures from 
the Hardy-Weinberg equilibrium was estimated using Y? test. 

Genetic distances "(D n and Lie, have been calculated accord- 
ing to both Ners (1972) and Rogers’ (WRIGHT 1978) methods. 
Calculations were performed using BIOSYS-1 computer program 
(SWOFFORD & SELANDER 1989). Times of evolutionary divergence 
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(t) were calculated using the formula proposed by NEI (1975) when 


(ELE 


t=5-10°D,.. 


A principal component analysis (PCA) was carried out using 
SYSTAT 5.2 (WILKINSON et alii 1992). The PCA was based on a 
matrix of the allele frequencies of polymorphic loci, after removing 
alleles with frequencies lower than 10%. 


Table 1 - Code number of Speleomantes strinatii samples listed with the collecting 


Number 


Soon Awe oS we 


es | 
— > 


— 
ww 


sites and their province and region. 


Collecting sites 


Pietra di Vasca Mt. 

Eastern Slope of Groppi Mt. 

Environs of Carro and Ponte S. Margherita 
Pù Mt., in the Rio Frascarese Valley 

Near Cassagna 

Cave “Tana de Strie” 


Cave “Tana da Cruxetta” 


Near Maissana and SW slope of Baralucco Mt. 


Near Codolo 

Near Bardi 

Near S. Pietro di Novella and Rapallo 
Near Scaglia, Bargagli 

Near Cartasegna 

Near Isoverde 

Near Millesimo 

Near Roveirola, N slope of Bric Cornavento 
Near Finale Ligure 

Near Gavone 

Near Toirano 

Left bank of the Rio Roburentello stream 
Near Roaschia and near Rivoera 

S. Bartolomeo Hill, S slope 

Near Caravonica 

Downstream T'enarda dam 

Near Rocchetta Nervina 

Near Tenda 

Near Luceram 

Between Peille and St. Martin de Peille 


Near Aspremont 


Province 


Genoa 
Genoa 

La Spezia 
Genoa 
Genoa 

La Spezia 
La Spezia 
La Spezia 
Massa Carrara 
Parma 
Genoa 
Genoa 
Alessandria 
Genoa 
Savona 
Savona 
Savona 
Savona 
Savona 
Cuneo 
Cuneo 
Imperia 
Imperia 
Imperia 


Imperia 


Region/Country 
Liguria 
Liguria 
Liguria 
Liguria 
Liguria 
Liguria 
Liguria 
Liguria 
Tuscany 
Emilia Romagna 
Liguria 
Liguria 
Piemonte 
Liguria 
Liguria 
Liguria 
Liguria 
Liguria 
Liguria 
Piedmont 
Piedmont 
Liguria 
Liguria 
Liguria 
Liguria 
France 
France 
France 


France 
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RESULTS 


Eleven loci out of the 33 studied were monomorphic, show- 
ing the same allele in all the populations sampled: [dh-2, NADH- 
dh, Np, Sod-1, Sod-2, Aat-2, Adk, Pep-C, Lap, Fum, Mpi. The 
remaining 22 loci showed from two to four alleles per population. 
The pattern of allele frequencies at many loci concurred in show- 
ing a higher genetic variation in western populations, while samples 
from the eastern part of the range were genetically more homo- 
geneous. This pattern is confirmed by the distribution of the private 
alleles (i.e., alleles found in a single population; NEEL 1973): 15 over 
18 private alleles where in the samples from the western part of the 


La Spezia E 


Factor 2 


Factor 1 


Fig. 1 - A) Geographical distribution of the studied populations, symbols refer 
to the subspecies to which they were ascribed on morphological ground. 
B) PCA based on allele frequencies, the percentage of the total variance 
explained is 88.7 (80.6 is explained by Factor 1). The terra typica of each 
subspecies is circled. 


m (squares) S. s. ligusticus, ® (stars) S. s. argentatus and “nameless” samples 
from Finalese (15-18), & (triangles) S. s. bonzanoi, @ (circles) S. s. strinatii. 
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Table 2 - Parameters of genetic variability in the 29 studied populations of S. stri- 
nati. Standard errors are in parentheses. 


Population A P P Hi H 


99 95 o e 
i 1.1 Dai 6.1 0.009 0.011 
(0.005) (0.006) 
Da 1i 9.1 9.1 0.047 0.035 
(0.029) (0.020) 
ii m 6.1 6.1 0.027 0.023 
(0.019) (0.016) 
4, 1.1 9.1 9.1 0.022 0.022 
(0.013) (0.013) 
5, sve 6.1 3.0 0.010 0.009 
(0.008) (0.007) 
6. 1.1 6.1 3.0 0.007 0.007 
(0.006) (0.005) 
da 1.1 6.1 6.1 0.018 0.022 
(0.012) (0.016) 
8. wi 10M lo 0.034 0.028 
(0.020) (0.015) 
9. 1.1 6.1 6.1 0.014 0.019 
(0.010) (0.015) 
10. imo 3.0 3.0 0.013 0.011 
(0.013) (0.011) 
11. da 9.1 9.1 0.031 0.037 
(0.018) (0.022) 
ea “xi 9.1 9.1 0.035 0.030 
(0.022) (0.018) 
13. ii 9.1 9.1 0.025 0.027 
(0.016) (0.017) 
14. 1.1 9.1 9.1 0.027 0.024 
(0.017) (0.015) 
15. iat 9.1 9.1 0.036 0.031 
(0.022) (0.018) 
16. T Da DARI 0.079 0.085 
(0.030) (0.031) 
17. 1.2 18.2 18.2 0.062 0.056 
(0.028) (0.024) 
18. 9) 52 132 0.052 0.050 
(0.025) (0.024) 
19. 14 5 18.2 0.079 0.077 
(0.028) (0.027) 
20. 1.1 12 Tari 0.041 0.043 
(0.022) (0.023) 
21. tel 6.1 6.1 0.020 0.021 


(0.015) (0.015) 
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Table 2 - (Continued) Parameters of genetic variability in the 29 studied popula- 
tions of S. strinattt. Standard errors are in parentheses. 


Population A IAA A RE H, 
221 1.4 36.4 30.3 0.093 0.1125 
(02027 (0.036) 
ZO rS 2785 27.3 0.108 0.108 
(0.034) (0.034) 
24. 1.4 30.3 22 0.088 0.091 
(0.030) (0.031) 
25, k3 24.2 24.2 0.093 0.096 
(0.031) (0.031) 
26. 152 24.2 24.2 0.091 0.085 
(0.032) (0.030). 
D. 12 18.2 18.2 0.081 0.075 
(0.032) (0.030) 
28. 183 212 1522 0.048 0.065 
(0.020) (0.026) 
29: is O73 24.2 0.077 0.089 
(0.026) (0.030) 


range (samples 15-29), while only three have been found in the 14 
easternmost populations. Moreover, all the private alleles with fre- 
quencies higher than 10% have been found in western populations. 
The higher genetic variation of westernmost populations is also 
evidenced by the values of genetic variability parameters (Tab. 2). 
In the easternmost populations (1-14) the averaged expected hetero- 
zygosity was 0.022 (range 0.007-0.037) and the percentage of poly- 
morphic loci (P,,) had a mean value of 8.0% and a range of 3-12.1%. 
The values observed for the samples 15-21 were higher, with mean 
values of H, reaching 0.052 (range 0.021-0.085) and P, between 
6.1-33.3 (mean value 16.5). The western populations (samples 22- 
29) show the highest levels of genetic variability, with heterozygosity 
values ranging between 0.065 and 0.125 and a percentage of poly- 
morphic loci from 18.2 to 36.4. Mean values of these parameters 
were, hp —U,092 aneer,.—20.1, two to four-tumes those recorded for 
the other populations. 

The PCA evidences three main groups of samples (fig. 1B). 
The group in the upper-right part of the graph clusters the east- 
ern samples, which are very similar to each other, in spite they are 
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scattered over a half of the whole range of the species. This group 
corresponds to the subspecies ,S. s. ligusticus (Stefani, 1969), being 
the terra typica of this taxon Valdettaro Cave, near Rapallo (sample 
11). A second group was represented by the westernmost popula- 
tions (samples 22-29), including the samples (22-23) so far described 
as S. s. bonzanot (Bruno & Bologna 1973) and S. s. strinatii (Aellen, 
1958) (samples 24-29). The samples belonging to the subspecies S. 
s. argentatus (Stefani, 1969) (samples 19-21) and those belonging to 
a nameless form (15-18) were split in two differentiated groups. The 
first group (samples 19-21) was close to S. s. ligusticus. The second 
group was in the lower part of the graph and includes the Ligu- 
rian populations from the “Finalese” area (samples 15-18). These 
four populations are the most genetically differentiated in spite they 
are located in the central part of the species range. Therefore, PCA 
indicates that three main population groups could be found within 
S. strinatit: western populations (22-29), samples from “Finalese” 


(15-18), and eastern populations (1-14 and 19-21). 


Table 3 - Genetic distances calculated according to Rogers’ method (WRIGHT, 
1978) for the three population groups identified within S. strinatii. Mean 
values below the diagonal, range of distances above the diagonal. On the 
diagonal are the ranges within each group. 


Population group “Eastern” “Finalese” “Western” 


(samples 1-14,19-21) 
(samples 15-18) 
(samples 22-29) 


The genetic distances calculated according to Rogers’ method 
(WRIGHT 1978) supported the identification of three population 
groups (Table 3). The eastern samples (1-14 and 19-21) are shown 
to have a quite low genetic differentiation, with the easternmost 
populations (samples 1-14) showing very low DRogers (between 
0.008 and 0.141), that slightly grew to 0.127-0.212 when the three 
samples from the central part of the range (19-21) are added. The 
populations from western range (samples 22-29) are well differenti- 
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ated both from the other population groups (DRogers between 0.196 
and 0.381) and among them (ranging 0.108-0.344), confirming the 
genetic heterogeneity of this cluster. Finally, the populations 15-18 
from “Finalese” are confirmed as highly differentiated, showing the 
highest recorded average distance from the other populations (aver- 
age DRogers 0.312). 


Table 4 - Mean genetic distances (Dyp NEI, 1972, below the diagonal) and diver- 
gence times in years (t, NEI, 1975, above the diagonal) between the three 
genetically characterised forms identified within S. strinatit. 


“Eastern” 
“Eastern” form (samples 1-14,19-21) t = 590.000 t = 375.000 


“Finalese” form (samples 15-18) De. : t = 500.000 


Nei 


“Western” form (samples 22-29) aaa 


DISCUSSION AND CONCLUSIONS 


The results obtained allowed the identification of three geneti- 
cally and geographically characterised groups of populations. A 
genetically homogeneous group occupies a wide area (samples 1-14 
from the eastern part of the range) and is genetically close to the 
population group from the western area (19-21). The relatively small 
westernmost part of the range hosts populations (22-29) that are 
highly differentiated both among themi and from the other clusters 
identified. Finally, the most divergent population group is located 
in the central part of the range and covers a very small area in the 
“Finalese” region (samples 15-18). 

From a taxonomic point of view, it was hard to find out a cor- 
respondence between these three groups and the formally described 
morphological subspecies. The “eastern” group (1-14 and 19-21) 
includes the terra typica of the subspecies S. s. ligusticus (Stefani, 
1969), but also that of S. s. argentatus (Stefani, 1969) (cave Grotta 
delle Ruccaie, near Toirano, sample 19), along with the populations 
from Roburent area, usually included in this last subspecies. In the 
westernmost group are clustered all together the populations classi- 
fied as S. s. strinati (Aellen, 1958) (with the terra typica: the cave 
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Grotte d’Aspremont, n. 29) and S. s. bonzanoi (Bruno & Bologna, 
1973) (samples 22 and 23, with this latter as the terra typica from 
the cave Tana I° du Casa). The most genetically differentiated clus- 
ter from the Finalese area does not include any terra typica. ‘These 
results clearly indicate that the so far proposed subspecies are incon- 
sistent with the genetic structure of the species. 


The genetic subdivision of S. strinati in the three above iden- 
tified population groups was consistent with the paleoclimatic his- 
tory and the orography of the species’ range. It is known that the 
speciation process in Speleomantes occurs mainly by divergence in 
allopatry, as in the larger part of plethodontids (GARC{iA-PARIS et 
alii 2000), and is strongly influenced by paleoclimatic events via the 
rising of geographic barriers allowing the differentiation via muta- 
tion sum (FORTI et alzi 1998; CIMMARUTA et alti 1998). In particular, 
it has been speculated that the glaciations played an important role 
at the evolutionary level, allowing the origination of the extant spe- 
cies by means of the alternation of glacial and interglacial periods 
during Quaternary (FORTI et alzi 1998). 


To test for the influence of glacial events on the genetic struc- 
turing of S. strinatit, the genetic distances between the so far identi- 
fied groups of populations have been used to estimate the divergence 
times among them, according to Nei’s calibration (NEI 1975). The 
obtained times have been compared with the paleoclimatic events 
occurred in the area, in particular with glacial/interglacial stages, in 
order to determinate if glacials can be identified as subdiving events, 
promoting population fragmentation and producing the extant ar- 
ticulated genetic structure*of the species. These results are reported 


in Tab. 4. 


The divergence of the most differentiated group from the Fina- 
lese area started during glacial stages: during Gtinz with respect to 
the “eastern” group (about 590.000 ya) and during Mindel (about 
500.000 ya) with respect to the “western” group. In the Mediter- 
ranean area, the glacial stages were characterised by a dry and cold 
climate, while the interglacials had a more mild and humid climate 
(WEST 1969; BERTOLANI-MARCHETTI 1985; PIGNATTI 1994). As a 
consequence, arboreal vegetation and animal species were restricted 
in warmer refugia at low altitude and along coastal areas during 
the glacials. The populations from the Finalese live in a small area, 
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completely surrounded by relatively high mountains (Mt. Settepani, 
Mt. Carmo, Mt. Beigua), that have probably acted as efficacious 
geographic barriers during ice ages, causing isolation and differen- 
tiation of the local remnant population group. 


The divergence of the two most far apart population groups 
(westernmost vs. easternmost) cannot be reconnected to a glacial 
stage. It is dated 350-400.000 years ago, at an interglacial stage, 
when higher altitude sites were colonised due to the range expan- 
sion of the arboreal vegetation. The genetic pattern of the “eastern” 
population group suggests a relatively recent range expansion, as 
witnessed by its high genetic homogeneity. Moreover, the close simi- 
larity with the populations from the Roburent area (samples 19-21), 
suggests an Apennine colonisation route from Emilia and Tuscany 
to western Liguria. The slightly higher genetic variability observed 
in the Roburent samples could be due to the “secondary contact” 
between the re-established populations invading the area and those 
resident in the Maritime Alps. 


This scenario is supported by the observed pattern of allele fre- 
quencies in the populations from Roburent area, holding “western 
alleles” together with the more common “eastern alleles”. Moreover, 
the Maritime Alps are a well known suture zone for either plant and 
animal species and the depicted pattern has been proposed for many 
other taxa (KROPF et alzi 2002 and references herein). This sce- 
nario could explain the impossibility of dating the divergence time 
between easternmost and westernmost population groups, being the 
data affected by the subsequent contact and gene flow between these 
two forms. 
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ABSTRACT 


The genetic structure of the plethodontid salamander Speleomantes strinatii 
(Aellen, 1958) has been analysed by means of allozyme electrophoresis. The studied 
samples were taken from the whole range of the species (north-western Italy and 
south-eastern France), to assess the degree and the pattern of genetic differentiation 
and to find out possible correlation with the paleoclimatic history of the area. The 
results obtained showed three genetically characterised population groups, whose 
distribution and levels of differentiation were related to the quaternary alternation of 
glacial and interglacial stages in the region. The subdivision of this species in four 
morphological subspecies is not supported by the genetic structure here evidenced. 
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RIASSUNTO 


Gli effetti delle glaciazioni quaternarie sulla struttura genetica di Speleo- 
mantes strinatii (Aellen, 1958) (Amphibia, Plethodontidae). 


La struttura genetica del geotritone Speleomantes strinatii (Aellen, 1958) è stata 
analizzata mediante analisi elettroforetica di sistemi gene-enzima, su popolazioni 
campionate in tutto l’areale della specie. Il grado di eterogeneità inter- ed intrapopo- 
lazionale, i livelli di variabilità e le relazioni di affinità genetica tra le varie popolazi- 
oni sono stati confrontati con la storia paleoclimatologica dell’area per configurare 
1 possibili processi microevolutivi che hanno originato l’attuale struttura genetica 
della specie. I risultati ottenuti hanno evidenziato la presenza di tre gruppi di popo- 
lazioni geneticamente ben caratterizzati, la cui origine è probabilmente riconduci- 
bile all’influenza dell’alternarsi dei periodi glaciali ed interglaciali del Quaternario. 
I dati ottenuti sono stati anche utilizzati per verificare la validità delle sottospecie 
morfologiche attualmente riconosciute, che non risultano essere sostenute da basi 
genetiche ed evolutive. 


